Effect of superradiance on transport of diffusing photons in cold atomic gases 
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We show that in atomic gases cooperative effects like superradiance and subradiance lead to a 
potential between two atoms that decays like 1/r. In the case of superradiance, this potential is 
attractive for close enough atoms and can be interpreted as a coherent mesoscopic effect. The 
contribution of superradiant pairs to multiple scattering properties of a dilute gas, such as photon 
elastic mean free path and group velocity, is significantly different from that of independent atoms. 
We discuss the conditions under which these effects may be observed and compare our results to 
recent experiments on photon transport in cold atomic gases. 
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The issue of coherent multiple scattering of photons in 
cold atomic gases is important since it presents a path 
towards the onset of Anderson localization transition, 
a long standing and still open issue. The large reso- 
nant scattering cross-section of photons reduces the elas- 
tic mean free path to values comparable to the photon 
wavelength for which the weak disorder approximation 
breaks down thus signaling the onset of Anderson lo- 
calization transition [1, 2]. Another advantage of cold 
atomic gases is that sources of decoherence and inelastic 
scattering such as Doppler broadening can be neglected. 
Moreover, propagation of photons in atomic gases differs 
from the case of electrons in disordered metals or of elec- 
tromagnetic waves in suspensions of classical scatterers 
for which mesoscopic effects and Anderson localization 
have been thoroughly investigated [1]. This problem is 
then of great interest since it may raise new issues in the 
Anderson problem such as change of universality class 
and therefore new critical behavior. New features dis- 
played by the photon-atom problem are the existence of 
internal degrees of freedom (Zeeman sublevels) and coop- 
erative effects such as subradiance or superradiance that 
lead to effective interactions between atoms [3]. These 
two differences may lead to qualitative changes of both 
mesoscopic quantities and Anderson localization. Some 
of the effects of a Zeeman degeneracy have been investi- 
gated in the weak disorder limit [4] using a set of finite 
phase coherence times [5] which reduce mesoscopic ef- 
fects, such as coherent backscattering [1, 6]. The aim of 
this paper is to investigate the influence of cooperative 
effects and more specifically of superradiance on the mul- 
tiple scattering of photons. We show that two atoms in 
a Dickc superradiant state [7] interact by means of a po- 
tential which, once averaged over disorder configurations, 
is attractive at short distances and decays like 1/r. This 
potential, analogous to the one considered in [8, 9], has 
important consequences on transport properties since the 
contribution of superradiant pairs of atoms in a dilute gas 
provides smaller values of both group velocity and diffu- 
sion coefficient so that the photons become closer to the 
edge of Anderson localization. 



Atoms are taken as degenerate two-level systems de- 
noted by \g) — \j g — 0, m g — 0} for the ground state and 
|e) = \j e = l,m e ) for the excited state, where j is the 
total angular momentum and m is its projection on the 
quantization axis, taken as the z axis. The energy sep- 
aration between the two levels including radiative shift 
is ftioo and the natural width of the excited level is HT. 
We consider a pair of such atoms in an external radiation 
field and the corresponding Hamiltonian is H = Hq + V, 
with 



Ho = n^o^(\e)(e\-\g)(g\) l + 



1=1 



(1) 



Q,k e (resp. a^ e ) is the annihilation (resp. creation) oper- 
ator of a mode of the field of wave vector k, polarization 
£k and angular frequency u)% = c|k|. The interaction V 
between the radiation field and the dipole moments of 
the atoms may be written as 



V = -dx ■ E(n) 



E(r 2 ) 



(2) 



where di is the electric dipole moment operator of the 
l-th atom and E(r) is the electric field operator. 

The absorption of a photon by a pair of atoms in 
their ground state, leads to a configuration where the 
two atoms, one excited and the second in its ground 
state, have multiple exchange of a photon, giving rise 
to an effective interaction potential and to a modi- 
fied lifetime as compared to independent atoms. These 
two quantities are obtained from the matrix elements 
of the evolution operator U(t) between states such as 
|<7ie2;0). There are six unperturbed and degenerate 
states with no photon, given by {|<7ie2i;0), |eij<72; 0)} in 
a standard basis where i,j = —1,0,1. The symmetries 
of the Hamiltonian, namely its invariance by rotation 
around the axis between the two atoms, and by reflec- 
tion with respect to a plane containing this axis, allows 
one to find combinations of these states that are given 
by \<t>\) = 75 [l e H32;0) + e|gie 2 ,;0)] with e = ±1, so that 
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(<f>* MtM) = 6 ij 5 ee ,St(t) and 

= (e lig2 ;0\U(t)\e lig2 ;0) + e(g 1 e 2l ;0\U(t)\e ll g 2 ;0) 

(3) 

The states \(f>\) are the well-known Dicke states, otherwise 
defined as \LM), where L is the cooperation number and 
M is half of the total atomic inversion [7] so that \4>f) — 
1 10) and \<j\) = 1 00) . For large times, t 3> r/c, where r is 
the distance between the two atoms, up to second order 
in the coupling to the radiation, we obtain that 



S«(t)~l-|[ABf-i^]. 



(4) 



The two real quantities AEf and T\ are respectively the 
interacting potential and the probability per unit time of 
emission of a photon by the two atoms in a Dicke state 
A standard calculation [10] gives 



AEf 



3HT 



COS fci 



or 



knr 



and 

r £ 3 
is. = l--c 
r 2 



-Pi- 



sin kor 
k a r 



cos kor sin kor 
Jko~rJ + JkoW 



sin kor cos kor 
JkoW ~ Jko~rJ 



(5) 



(6) 



where fco = i^o/ c - We have defined pi = 1 — rf and 
qi = I — 3??, f being a unit vector along the two atoms. 
At short distance fcor <C 1, we obtain that = 2Y for 
the superradiant state = 1 10) and = for the 
subradiant state = 1 00) . 

For a photon of wavevector k incident on an atomic 
cloud, the potential we shall denote by V e is obtained by 
averaging upon the random orientations of the pairs of 
atoms. Since (qi) = and (pi) = 2/3 regardless of i, we 
obtain for the average potential V e 



eV e (r) = (AEt) = -e 



hT cos kor 
~~2 k a r 



and the average inverse lifetimes of Dicke states are 



l + e 



sin kor 
k r 



(7) 



(8) 



which retains the same features as (6) for k r <C 1. 

Let us characterize the interaction potential V e . 
Whereas for a single pair of atoms, the potential (5) 
is anisotropic and decays at short distance like 1/r 3 , a 
behavior that originates from the transverse part of the 
photon propagator, we obtain that on average over angu- 
lar configurations, the potential (7) between two atoms 
in a Dicke state M = becomes isotropic and decays like 
1/r. This behavior is also obtained by considering the 
interaction of two-level atoms with a scalar wave. This 
could have been anticipated since the transverse contri- 
bution qi to the photon propagator averages to 0. A 
similar expression for the interacting potential has been 



obtained for the case of an intense radiation field [8, 9]. 
But this latter potential is fourth order in the coupling to 
the radiation and it corresponds to the interaction energy 
between two atoms in their ground state in the presence 
of at least one photon. The average potential V e we have 
obtained is different. It is second order in the coupling to 
the radiation and it corresponds to the interaction energy 
of Dicke states M = in vacuum. 

We turn now to scattering properties of Dicke states. 
The collision operator is given by T(z) = V + VG(z)V 
where V is given by (2) and G(z) is the resolvent whose 
expectation value in the Dicke state M = is ob- 
tained by a summation of the series of exchange of a 
virtual photon between the two atoms. The matrix 
element that describes the transition from the initial 
state |i) = |1 — l;ke) where the two atoms are in their 
ground state in the presence of a photon (ke) to the fi- 
nal state [/) = |1 — l;k e ) is the sum of the superradi- 
ant and subradiant contributions, T = T + + T~, with 
T± = (f\V\^)(^\G(u -loo^^Mi) [11]. A stan- 
dard derivation leads to the following expressions for the 
average amplitudes Tf- 



T + _ ^ e j(k-k')-R 



cos 



k • r 



cos 



G+ (9) 



and 



27 = Ae«( k - k ') R sin (^-) sin ( K^) G~ . (10) 



Wc have defined r = ri — r 2 , R = (ri + r 2 )/2 and 
A = j^d 2 (ej ■ e)(e* ■ i'*) (d is a reduced matrix element 
and Q the quantization volume). The average propaga- 
tors Gf associated respectively to the superradiant and 
subradiant states are, 



Gt ^ (^\G{5)\^) = 



(ii) 



where close to resonance, 5 = u) — uoo *C and where 
we have used (7) and (8) for the average potential and 
for the average inverse lifetimes. At short distances 
kor <C 1, the subradiant amplitude T~ becomes negli- 
gible as compared to the superradiant term (9). There- 
fore, the potential (7) is attractive and decays like 1/r. 
We can interpret these results by saying that, at short 
distances (kor <C 1), the time evolution of the initial 
state IV'(O)) = |ei,g 2 ;0) = ^7j[|0 + ) + I0~)] corresponds 
for times shorter than 1/r to a periodic exchange of a 
virtual photon between the two atoms at the Rabi fre- 
quency ((AE~) - (AE + ))/h ~ r/(fc r) which is much 
larger than T. For larger times, the two atoms return to 
their ground state and a real photon (k'e') is emitted. At 
large distances (kor » 1), the Rabi frequency becomes 
smaller than T, so that the excitation energy makes only 
a few oscillations between the two atoms, thus leading to 
a negligible interaction potential [12]. 
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It is interesting to derive the previous results in another 
way that emphasizes the analogy with weak localization 
corrections [1, 2]. 




FIG. 1: Diagrammatic representation of the two amplitudes 
that describe double scattering of a scalar wave. The wavy 
line accounts for the photon exchange between the two atoms. 
This diagram is known in quantum mesoscopic physics as a 
Cooperon. 

To that purpose, we write the scattering amplitude T 
defined previously as a superposition of two scalar ampli- 
tudes T\ and T 2 [13], each of them being a sum of single 
scattering and double scattering contributions, that is 



Ti = 



t 2 G 2 



i(k-k')-n 



ta 



i(k-ri— k'-r 2 ] 



and 



1 



t 2 G 2 



e i(k-k')T 2 +tGoe *(k.r 2 -k'. ri ) 



(12) 



(13) 



Here t = (2irT / ho) / (8 + iT / 2) is the amplitude of a scalar 
wave scattered by a single atom at the origin and the pref- 
actor t/{\ - t 2 Gl) where G = ~e tk ° r /inr accounts for 
the summation of the series of virtual photon exchange 
between the two scatterers. We single out in the total 
amplitude T = T\ + T 2 the single scattering contribu- 
tion T s and write the intensity associated to the double 
scattering term shown in Figure 1 as 



\T-TJ 2 = 2 



t 2 G Q 



1 



t 2 G\ 



[1 + cos(k + k') • (n - r 2 )] 



(14) 

We recognize in the bracket the well-known Cooperon in- 
terference term which is at the basis of coherent effects in 
quantum mesoscopic systems such as weak localization 
and coherent backscattering [1, 2, 6]. The interference 
term reaches its maximum value 1 for ri = r 2 so that we 
obtain from (12) and (13) that Ti = T 2 oc (1/2)T+, up 
to a proportionality factor [13]. Thus, the total ampli- 
tude is exactly given by the superradiant term with no 
subradiant contribution. 

We consider now multiple scattering of a photon by 
superradiant pairs built out of atoms separated by a dis- 
tance r and coupled by the attractive interaction po- 
tential V e . This situation corresponds to a dilute gas 

— 1/3 

that fulfills r -C A <C n i where rn is the density 



of pairs and Ao = 2ix jko is the atomic transition wave- 
length. Based on this inequality, we may consider the two 
atoms that form a superradiant pair through exchange of 
a virtual photon as an effective scatterer and neglect co- 
operative interactions between otherwise well-separated 
pairs. The photon behavior is described by the configu- 
ration average of its Green's function whose expression 
is obtained from a standard derivation [1]. In the limit 
of large enough densities of weakly scattering pairs, it 
reduces to the calculation of a self-energy given in terms 
of the average propagator (11) by 



6irhrni —+ 



6irrii 



rrij f Tm dr 
r m Jo y+ 2^7 



+ i 



(15) 



The average, denoted by 7TT , is taken over distances r 
up to a maximal value r m -C fc^ 1 which accounts for all 
possible mechanisms that may break those pairs. In the 
considered limit, the density of the gas can be assimilated 
to that of the pairs. The imaginary part of defines 
the elastic mean free path l e by ko/l e = -Ims' 1 ' , namely 



1 

z 



37T7J,- 1 



k 2 



k r r , 



p2k r„ 
JO 



dx 



(16) 



It is interesting to compare l e to the mean free path 
lo = ^1 + (2S/T) 2 ^j that corresponds to near reso- 
nant elastic scattering of a photon by independent atoms. 
At resonance (S = 0), we have l f- = |(fcor m ) 2 <C 1. Away 



from resonance, the elastic mean free path l e becomes 
smaller than lo and for blue detuning it is reduced in a 
ratio roughly given by \/{kor m ) 2 . 

Another important physical quantity is the group ve- 
locity v g given in terms of the refraction index r\ by 



c/vg = 7] + WgJ. Since is proportional to the po- 
larizability, the refraction index depends on its real part, 

namely r\ = (1 — (d/a^RcE^) 1 / 2 . From (15), we notice 
that rj ~ 1 for all values of the detuning S/T and in a 
large range of densities m so that 

- = 1 --^— f{kor m ,L (17) 
v„ n c 2k r m T 



where we have defined n c 



6-7T LUq 



and the function 



f(k r m ,A) 



1\2 



2k r m 1 - (A , 

dx ^~ 

(1 + (A + I) 2 



2 • 



(18) 



This expression of v g diverges at a large and negative 
value of the detuning £ ~ — 2k ^ r and beyond it takes 
both positive and negative values. Otherwise it is well 
behaved, meaning that it remains finite and positive for 
all values of the density rij. At resonance, the group 
velocity is 

AlTTli LOo 



v„ 1+ B r 



{k r m ) . 



(19) 
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The present expression of v g differs substantially from 
the one obtained for light interaction with independent 
two-level atoms. There, for densities rij > n c where n c 
defined above is usually overwhelmingly small, the group 
velocity is known to diverge at two symmetric values of 
the detuning of order unity and takes negative values in 
between (i.e., also at resonance). For instance, in a gas 
of Rb 85 atoms, where n, = 6 x 10 10 cto~ 3 , Ao = 780nm 
and ^ = 5.9MHz, we have ^ ~ 10 5 . The validity of 
the concept of group velocity in such systems has thus 
been often questioned [14] and an energy velocity has 
been defined which describes energy transport through a 
diffusive medium [15]. 

Transport of photons through a diffusing gas is charac- 
terized by the diffusion coefficient D = \v g l e that com- 
bines the elastic mean free path and the group velocity 
[1, 16], both derived from the complex valued self-energy 
(15). The diffusion coefficient D is of great importance 
since it enters in expressions of all measured physical 
quantities such as reflection and transmission coefficients, 
angular correlations of speckle patterns, time correlation 
functions of the intensity (diffusing wave spectroscopy), 
etc. [1]. Moreover, the critical behavior of transport close 
to Anderson localization transition at strong disorder is 
also obtained from the scaling form of D. Its expression, 
deduced from (16) and (17), depends on the range r m 
and on the detuning 6/T. Since the group velocity and 
the elastic mean free path are significantly modified for 
the case of superradiant pairs, we thus expect the diffu- 
sion coefficient to be different from its value obtained for 
independent atoms. We also define the transport time by 
i~tr(d) = le/vg = ^r- At resonance, it can be rewritten 

with the help of (19) as r tr (0) = which is consistent 
with our considering of superradiant pairs. We would 
like nevertheless to call attention to the fact that, away 
from resonance, Tt r (S) depends on frequency. 

We now compare our results to recent measurements 
of the diffusion coefficient D and of the group velocity v g 
obtained for multiple scattering of light at resonance, in 
a cold atomic gas of Rb 85 [17]. Since the range r m cannot 
be directly determined, we first use Eqs. (16) and (17) to 
obtain an expression independent of k r m given by the 
ratio %^ = ^ = 2r/c 2 . For Rb 85 atoms, this ra- 
tio equals 8.2 x 10~ 10 s/m 2 , which is in good agreement 
with the value 4.8 x lCT 10 s/m obtained from measure- 
ments of D and v g . Finally, from the previous numerical 
expression we deduce for the maximal range of interac- 
tion r m the value kor m ~ 0.51 also consistent with our 
assumption of superradiant states. Therefore, multiple 
scattering of photons by superradiant pairs provides a 
relevant mechanism that needs to be considered, in ad- 
dition to others e.g. scattering by independent atoms, 
for description of multiple scattering properties of dilute 
cold atomic gases. 

We have considered multiple scattering of a photon on 



pairs of atoms that are in a superradiant state. On aver- 
age over disorder configurations, an attractive interaction 
potential builds up between close enough atoms that de- 
cays like 1/r. The contribution of superradiant pairs, 
resulting from this potential, to scattering properties is 
significantly different from that of independent atoms. It 
leads to a well defined but much smaller group velocity as 
compared to c and correlatively to a smaller diffusion co- 
efficient. For densities considered in recent experiments 
on cold Rb 85 , the quantity kol e that describes eventually 
the closeness to a localization transition, is reduced at 
moderate detunings, by one order of magnitude. This 
effect is expected to be even stronger for larger densities 
which could then be close to the localization edge. This 
research is supported in part by the Israel Academy of 
Sciences and by the Fund for Promotion of Research at 
the Tcchnion. 
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